Ramsey interference in one dimensional systems: The full distribution function of 
fringe contrast as a probe of many-body dynamics 
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We theoretically analyze Ramsey interference experiments in one dimensional quasi-condensates 
and obtain explicit expressions for the time evolution of full distribution functions of fringe contrast. 
We show that distribution functions contain unique signatures of the many-body mechanism of 
decoherence. We argue that Ramsey interference experiments provide a powerful tool for analyzing 
strongly correlated nature of ID interacting systems. 



Introduction Recent progress in the field of ultracold 
atoms not only expanded our understanding of equilib- 
rium properties of interacting Id Bose gases [J, |3| but 
posed new theoretical challenges by studying far-from- 
equilbrium dynamics of such systems. Recent experi- 
ments addressed such questions as thermalization and 
integrability Isj , decoherence after the splitting of two 
condensates [j] and spin dynamics of two component Bose 
mixtures [1]. Motivation for such experiments comes not 
only from the basic interests in many-body dynamics [6^| 
but also from possible applications of ultracold atoms 
such as quantum information processing j?! and intcr- 
ferometric sensing J^. In this paper we theoretically an- 
alyze the decoherence dynamics of Ramsey interference 
fringes in one dimensional quasi-condensates. Such sys- 
tems have been considered for possible applications in 
atomic clocks and quantum enhanced metrology 0, [lo[ . 
In this paper we show that Ramsey interferometer is also 
a powerful tool for studying many-body dynamics of low 
dimensional quantum systems. We find that decoherence 
of Ramsey fringes is strongly affected by the multimode 
character of one dimensional systems. Moreover we will 
demonstrate that time evolution of the full distribution 
function of fringe contrast provides unique signatures of 
this many-body decoherence mechanism flTl. The idea 
of using noise and distribution functions to character- 
ize equilibrium many-body states of ultracold atoms has 
been discussed in several theoretical papers and ap- 
plied in experiments [l|, However there has been so 
far no application of this approach to non-equilibrium 
dynamics. This paper constitutes the first proposal to 
study non-equilibrium dynamics of ultracold atoms with 
quantum noise. 

The role of interactions in Ramsey interferometers 
with BEG was first addressed in the pioneering paper 
of Kitagawa and UedaH. They used single mode ap- 
proximation to predict the interaction induced deco- 
herence of Ramsey fringes along with the appearance 
of spin squeezed states. Their work stimulated ideas 



for quantum-enhanced metrology that take advanta ge o f 
spin squeezed states formed in interacting BECs (ol Il4|. 
For the analysis of one dimensional quasi-condensates, 
however, the single mode approximation cannot be ap- 
plied because these systems do not have macroscopic oc- 
cupation of a single state even at zero temperature. The 
non mean-field character of the multi-mode spin dynam- 
ics in ID quasi-condensates was first reported in the ex- 
periments of Widera et al However, this work did 
not provide the definitive demonstration of the many- 
body origin of decay. In the following, we argue that 
unambiguous signatures of the multimode decoherence 
are found in the full distribution function of the Ramsey 
fringe amplitudes. Such distribution functions should be 
accessible in experiments with ID quasi-condensates re- 
alized on Atom Chips ^iB'l, because such systems do not 
average over multiple tubes and thus allow the measure- 
ments of shot-to-shot fluctuations [if. 

Now we describe the Ramsey sequence considered in 
this paper. Here we identify two hyperfine states as spin 
up and down states. Ramsey sequence is carried out as 
follows: (i) all spins of the atoms are prepared in the spin 
up state; (ii) 7r/2 pulse is applied to rotate each spin into 
the X direction; (iii) spins freely evolve (precess) for time 
i; (iv) another ■n/2 pulse is applied to map the trans- 
verse spin component to the z direction, which is then 
measured. Measurements yield a net spin Si for a seg- 
ment of length I and we assume that / is smaller than the 
the system size but large enough to contain large number 
of particles iV; >> 1. In such a case, the simultaneous 
measurements of Si and 5*^^ are possible, because even 
though operators Sf and S'f generally do not commute, 
non-commutativity gives only corrections of the order of 
1/\/Ni relative to the average values Commuta- 
tivity of Si and Sf implies, in particular, that we can 
define the joint distribution function P^'^ for the two 
transverse spin components Sfand Sf . In experiments, 
measurements of Pf^'^ is possible by mapping the spin 
orientations in x — y plane to z direction by tt/2 pulse. 
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FIG. 1: (a),(c): Evolution of joint FDF Pj^'^ with short 
integration length Z/^s — 10(top, (a)) and long integration 
length l/^s = 30(bottom, (c)). (b),(d): Corresponding FDF 
for spin x, with short integration length //^s = 10 (top, 
b)) and long integration length l/£^s ~ 10(bottom, d)). Here 
L/£,s ~ 200, Ks — 20. Time is measured in units of Cs/cs. 



followed by local measurements of Sz [l3|- Sf and Sf 
as well as the magnitude of spin = a/ {S^)^ + {Sf)^ 
can be found by taking the integration over /. The an- 
alytic solution for the time evolution of Pj^'^ constitutes 
the main result of this paper. In addition, we assume 
that I is larger than the spin healing length , so we can 
use Tomonaga-Luttinger liquid approach to describe the 
collective spin dynamics (see also below). For simplicity 
we work in the rotating frame of the Larmor precession, 
and consider the spins before the last tt/2 pulse. Then 
the amplitude of Ramsey fringes, as it is conventionally 
defined, corresponds to S[. 

Our main results are summarized in Fig. [T] and can 
be understood from the following physical arguments. 
Strong fluctuations present in ID systems forbid the ex- 
istence of long range coherence [rs"], and spatial fluctua- 
tions coming from different wavelength strongly affect the 
dynamics in ID. Among those, fluctuations with wave- 
lengths longer than the integration length I rotate Si as 
a whole. So they decrease but not the magnitude 
of the spin ^ ^WF+WF- Fluctuations with 
wavelengths shorter than I decrease both 5*;^ and si- 
multaneously. Fig[T] (a) (b) show the situation where 
fluctuations with wavelength larger than I dominate the 
dynamics. In (a), we see that the magnitude of the spin 
decays only slightly from the initial state but the 
direction of the spin is randomized during the time evo- 
lution. Note that in this case the distribution function 
of 5f has a very peculiar shape with two peaks at large 
positive and negative values. We call this regime, "spin 
diffusion" regime. Fig[T] (c) (d) show the situation where 
fluctuations with wavelengths shorter than I dominate. 



In this case Sj^ and Sf decay in the same timescale. We 
call this regime, "spin decay" regime. Below we argue 
that the crucial parameter of the system is a dimension- 
less ratio proportional to the length of the integration 

2 T 

region = 4^ ^ • When Iq < 1 the system is in the spin 
diffusion regime and the other limit >> 1 is the spin 
decay regime. 

Model. Following the flrst 7r/2 pulse we have a two 
component Bose mixture with equal densities of both 
species. Tomonaga-Luttinger liquid (TLL) approach, 
which we use in this paper, focuses on the linearly dis- 
persing modes in the low energy part of the spectrum. 
For simplicity we consider the case when interaction pa- 
rameters satisfy = gi^i 21|. This condition can be 
reached for the hyperfine states \F = l,nip — —1) and 
\F = 2,tof = +1) of ^^Rb that are commonly used in 
experiments [H When this is the case, the charge 
and spin parts of the Tomonaga-Luttinger Hamiltonian 
decouple and the spin part of the Hamiltonian is given 
by 
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Here L is the total system size, ns{r) describes the lo- 
cal spin imbalance (i.e. z component of the spin) ng = 
^a{^'^ap)i^P and Ug^k is the Fourier transform of ns{r). 
(ps {r, t) describes the direction of the transverse spin com- 
ponent pe**^" = V'a(^ "'0/3)^/3 'with p being the average 
density for each species. Variables Us and <j)s obey canon- 
ical commutation relations [ns(r), 0s(^')] = —i5{r — r'). 
Ks is spin Luttinger parameter representing the strength 
of interactions 0, and Cs is spin- wave velocity. Hamilto- 
nian ^ has momentum cutoff set by the inverse of the 
spin healing length ^7^. In the weak interaction limit, 
these parameters are related to physical parameters as 
Cs = yjgsplm, 6 = and Kg ^ ^spf^, where m is the 



4tt h +°ij-~2°Tj- 



is the 



mass of the particle and gs — ^ 2 
interaction strength in spin channel. In this paper, we 
focus on the regime gs > Q when the system is miscible. 
Our approach can be extended to < but will be lim- 
ited to times before z magnetization per atom becomes 
of the order of one. The first and second term of ^ 
correspond to A: ^ and fc = part of the Hamiltonian, 
respectively. Here operators 6^ ^ create spin excitations 
with momentum fc, and these spin excitations are the 
main focus of our study. 

Transverse part of the spin operator Si is given by 

1-1/2 

,(r)), 5f = / drpsin(0,(r)),(3) 
J -1/2 

When describing spin dynamics one typically considers 
time evolution of the expectation values {S^{t)). How- 
ever important information is also contained in the shot 
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to shot fluctuations of S^{t). Such quantum noise is 
captured by full distribution functions (FDF) of spin 
operators,P;°(a, t) In particular, high moments of 

Sf{t) can be obtained from these FDFs Pl^{a,t). Physi- 
cally Pj^ia, t)da is the probability that a single measure- 
ment of the spin operator at time t gives the value 
between a and a + da. In the experiments, Pl^{a, t) can 
be obtained by making histograms of the measurement 
results oiSfit). 

To describe time evolution of spin operators ([3]) we 
need to characterize the initial state of the system after 
the first 7r/2 rotation. The difficult part is translating 
the initial state of the microscopic language, where spins 
of all atoms pointing in the x direction, to the one in 
terms of the coarse grained degrees of freedom (j^sif) and 
ns{r). Classically one expects the initial state to be the 
eigenstate of (j)s{T) with eigenvalue zero for all r. How- 
ever such state is unphysical in quantum mechanics since 
it leads to infinite uncertainty of the conjugate variable 
ns{r), and thus, to infinite energy. More sensible initial 
state is a squeezed state of harmonic oscillator ([T]) , which 
has reduced uncertainty in (ps at the expense of enhanced 
fluctuations in nc[5l.l20l| . To determine parameters of this 
state we observe that spins of individual atoms are inde- 
pendently rotated by the first 7r/2 pulse into x direction, 
so the initial state satisfies {S^{r)S^{r')) — ^5{r — r'). 
Thus we find a Gaussian state for the spin operator S"^ 
in momentum space with fluctuations p/2 for all k: 
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where 2Wk = (1 - Q:fe)/(1 -t- a^), ak = \k\Ks/Trp and TV 
is the overall normalization of the state. For the uni- 
form part of the spin operator we also have a squeezed 
state {nsfi\tps,k=o) = exp(— l/(2p)n^ q). We note that 
model ([T]) has a short distance cut-off so the S function 
in {S^{r)S^{r')) should be understood as rounded off on 
the scale of , which is implicit in the momentum cut-off 
in Eq. g]). 

Time evolution of the state (|4]) leads to Wk 
yY^^2iCs\k\t ^ From the resulting expression for the state 
at time t, one can readily calculate the decay of Ram- 
sey fringes given by {Sf{t)), which is independent of in- 
tegration length / (See also 0, [l^). To calculate time 
evolution of FDF, we deflne instantaneous annihilation 
operators ^ks{t) such that application of ^ks(t) on the 

state exp |^VFfce^*^=l'^l*5^ |0) gives zero. Using op- 

erators Tfcsft), one can apply the approach described in 
Refs. [UllH for calculating distribution functions of equi- 



librium systems. After direct calculation we find [21 1 
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where the real and imaginary part of a corresponds 
to X and y component of 5;, respectively. Eq. (O, 
Q allow a simple physical interpretation. Function 
x(r, t, {A}) defines the local direction of transverse mag- 
netization, which results from the summation over spin- 
wave like modes sin(fcr-|- Ae^fc). Amplitudes of individual 
modes are given by the time dependent expectation val- 
ues {\(j)s^k{t)\'^) and by the set of random variables Xrsk 
drawn from a Gaussian ensemble. Eq. ([5]), ([6|) reflect the 
key feature of dynamics of the quadratic Luttinger model 
(1): initial Gaussian state for (j)s,k remains Gaussian at 
all times [19J . 

Time evolution of (|(/)s_fc(t)p) following the flrst tt/2 
rotation can be understood as free dynamics of a har- 
monic oscillator. From the conjugate nature of (ps^k 

at 



and n^.fe we find (|0s,fc(O)r) = iTlTTTrMpy = Tp 
t — Q. Subsequently (|0s,fe(i)P) oscillates between the 
minimal value in the initial state and some maximum 
value {\(l)s,k\^)max at the frequency of a harmonic os- 
cillator Cs\k\. {\(t>s,k\'^)max Can bc estimated from en- 
ergy conservation. Since the initial state was squeezed 
with respect to 4>sk, most of the energy of the mode 
is stored in the interaction term |ns,fep. Therefore the 
total energy of the harmonic oscillator for momentum 
k can be approximated by which in turn gives 

{\(ps,k?)max ^ = ^ {^\^K:) ■ These considera- 

tions lead to the dynamics of phase fluctuation ampli- 
tude of the form in ([5]). We note that the spin fluc- 
tuations are dominated by small momentum since the 
maximum fluctuation amplitude is suppressed as l/Zc^ 
for large momentum. This justifles our analysis based on 
the Tomonaga-Luttinger theory. 

Results of numerical plots based on Eq. ^ are 
shown in Fig. [TJ In Fig. [2] (a) we also present the dis- 
tribution function Pj^ of the magnitude square of the 

integrated spin (S'^'^)^, which clearly demonstrates the 
difference between the " spin diffusion" and " spin decay" 
regimes. The character of these distribution functions 
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FIG. 2: (a):Time evolution of the distribution for the 
magnitude of spin (5*;^) • (b): time evolution of (Sf) and 

\l (('S'l^)^) with various integration length = 10,30,50. 
Here we set Ks=2Q, L/^s = 200. 



can be understood from the following arguments. We 
first discuss the "spin diffusion" regime, where the char- 
acteristic wavelength of spin fluctuations is longer than 
the integration length I (Fig. [T](a), (b)). In this regime, 
all spins within I essentially point in the same direction 
and Sj^ remains large even after a long time evolution. 

Thus we find a peak at {S/~)'^ « 1 in the distribution 
function P;^(red line). In the other regime of "spin de- 
cay" (Fig. [l](c), (d)), the typical length scale of spin fluc- 
tuations is shorter than the integration length I. In this 
case integration of spins over I is akin to taking a random 
walk in 2D plane and accordingly, distribution for {Sj^^ 

approaches exponential form with a peak at {S^)^ — 
in P;^(bluc line). In the intermediate regime, we observe 
two peak structure for Pj^, where the distribution ex- 
hibits both characteristic peaks (green line). We can un- 
derstand the condition that separates these " spin decay" 
and " spin diffusion" type dynamics in the following way. 
Deviation of spin angles aX r = I relative to r = can be 

: sin(fcZ). 

where 



jk ^rsk \/ {\^s,k \^)m,ax ' 



estimated from Ax w J2k 
A typical magnitude of Ax is given by ((Ax)^) 
the average is taken over fluctuations of Xrsk- The fac- 
tor sin(fc^) in Ax effectively limits momentum integra- 
tion range to A; > 27t/1, so ((Ax)^) ~ 2K g • When 
((Ax)^)^^^ is smaller than 2tt the system is in the "spin 
diffusion" regime. When ((Ax)^)""^^^ becomes of the or- 
der of 27T and larger, the system enters the "spin decay" 
regime. The crossover takes place around ^ ~ 1. 

In the "spin diffusion" regime, the dynamics of Sj^ and 
S'f display different time scales as can be seen in Fig. [5] 
(b) . In order to understand this separation of time scale, 
we note that the magnitude of the integrated spin Sj^ 
is only affected by fluctuations with short wavelengths, 
A < /, for which dynamics takes place at short time scale. 
Hence (5';-'-) decays until the time t± w 27r//cs(see Eq. 
([5])) and then it reaches a saturated value. On the other 



hand Sf is affected by excitations of all wavelength, so 
(Sf ) decays until it reaches 0. These behaviors are shown 
in Fig. [21 (b), where we compared the decay of (Sj^) for 
various segment length I and {Sf). 

Non-trivial time evolution of the distribution func- 
tions Pj^, P^'^ and especially the striking contrast of the 
"spin diffusion" and "spin decay" regimes should pro- 
vide unique signatures of the non-mean-field character 
and multimode dynamics of ID systems. 

Before concluding this paper we point out that our 
analysis can be extended to the problem of splitting a sin- 
gle ID quasi-condensate into two, as it was done with an 
RF potential in experiments reported in Refs.Q]. While 
earlier theoretical work focused on the time decay of the 
average fringe contrast [2^, the method developed in this 
paper can be used to study the time evolution of the full 
distribution function. 

Summary. We provided theoretical analysis of Ram- 
sey interference experiments with one dimensional quasi- 
condensates. We discussed time evolution of the full dis- 
tribution functions of fringe contrast and showed that 
they contain unique signatures of the many-body dynam- 
ics of one dimensional systems. This work was supported 
by the NSF grant DMR-0705472, Harvard MIT CUA, 
DARPA OLE program, AFOSR MURI, and Swiss NSF. 
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